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Abstract 

The string susceptibihty exponents of dynamically triangulated 2-dimensional 
surfaces with various topologies, such as a sphere, torus and double-torus, 
were calculated by the grand-canonical Monte Carlo method. These simula- 
tions were made for surfaces coupled to d-Ising spins ((i=0,l,2,3,5). In each 
simulation the area of surface was constrained to within 1000 to 3000 of trian- 
gles, while maintaining the detailed-balance condition. The numerical results 
show excellent agreement with theoretical predictions as long as d < 2. 
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1 Introduction 



In 2-dimensional quantum gravity one of the most fundamental physical (universal) 
quantities is the string susceptibility exponent (7) which characterizes the partition 
function of surfaces with a fixed topology through 

Z{A) ^ A^~3gA,A 

for the large-A limit, where A is the 2-dimensional volume (area), and Ac is the 
cosmological constant. The string susceptibility exponent is known to be given by 
the formula 

+ 2 (2) 

in the continuum framework [|l| as well as the matrix model [^, where h is the 
number of handles of the 2-dimensional surface, and c represents the central charge 
of matter fields coupled to the surface. Eq. (^ is expected to hold for a discretized 
2-dimensional surface with arbitral h and c < 1, when the dynamical triangulation 
recipe is adopted to generate random surfaces. 

In order to measure the string susceptibility exponent with any topology, we 
introduced the grand-canonical Monte Carlo method. Among quantities of interest 
the string susceptibility exponent of a spherical surface has been extensively calcu- 
lated numerically by various methods for examining the theoretical prediction, and a 
numerical method as well. The MINBU (minimum neck baby universe) algorithm 
and the grand-canonical Monte Carlo methods [^] have been successful as long as 
the central charge is less than one with sphere topology. However, in the MINBU 
algorithm measuring 7 for a surface with higher genus is more complicated than for 
a sphere, because the minimum neck baby universe must not be topologically non- 
trivial, and the baby universe distribution describing two topologically different 7 
must be separated independently 0. Also, this algorithm is subject to rather large 
finite-size effects, because the baby universes are dominated by small sizes, even if 
the mother universe is large. It cannot give reliable values for c 7^ unless we know 
the finite-size correction accurately 0. On the other hand, in the grand-canonical 
ensemble methods the finite-size effect enters through the total area, and becomes 
controllable in reasonably large size simulations. In this study we generated a suf- 
ficiently large size surface with a wider region of area compared with the previous 
studies, which is essential to determine Ac accurately. 

Eq.(^ shows that 7 becomes complex for 1 < c < 25 {h ^ 1), which is clearly 
unphysical and suggests that the nature of the surface will change significantly as 
c moves beyond unity. It is generally believed that this "transition" reflects the 
tachyonic instability of a bosonic string for c > 1. In a numerical simulation by 
dynamical triangulation the partition function is kept real for any value of c, and 
the instabihty is considered to appear along with the collapse of the surface to 
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branched polymers. The string susceptibihty exponent of the branched polymer 
ilBp) depends on the model. In fact, in ref. ■jep for a sphere is evaluated as 
|, while in ref. they calculated 'Jbp with any topology: 



13, 

7bp = - + -h. (3) 

In the case of a torus [h = 1) the effect of the tachyon cannot be seen in the string 
susceptibility exponent, because both eq.(@) and eq.(|) give the same value of 7 = 2 
for any c. 

In order to check the theoretical expectations (eqs.(||) and (||)), we need to 
calculate the string susceptibility exponent for a surface with any topology. 

We describe details of our method in Section ^, and our numerical results in 
Section |^. Section | is devoted to a summary and discussions. 

2 The method 

2.1 How to evaluate 7 

The grand-canonical partition function with a fixed topology is defined by 

Z{X)= Y.W{T,N2)e-''''^Z^, (4) 

T,N2 

where 

Z. = exp(/?f: ^afa^j. (5) 
[ /i=i <ij> ) 

In eqs.(D and (|^) T refers to a specific triangulation, W{T,N2) is the symmetry 
factor, Ni is the number of i-simplices {i.e. No, Ni and N2 are the number of 
vertices, links and triangles), A is a parameter, and /5 is a coupling strength of Ising 
spins. There are d kinds of Ising spins (o"f), and each Ising spin lives on a vertex 
of triangulated surface (not on the dual lattice). One kind of Ising spins carries 
i central charge at the critical point (3 = jSc, and thus the d kinds of Ising spins 
correspond to c = |. 

The partition function is written for the large- iV2 limit as 

Z(A) f^iJ^A^rV^^-^^^^^ (6) 
N2 

Then, the probability of fixing a surface with N2 triangles {P{N2)) for large N2 is 
approximately given by 

P(Ar2) oc iVr^e-(^-^^)^^ (7) 
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From 

if we tune the parameter A close to Ac within the precision 

|A - A, « (9) 

we can evaluate the string susceptibility exponent by plotting lnP(A^2) as a function 
of lnAr2. 

The parameter A is determined in the following manner. If A is chosen to be 
less (or greater) than Ac, the system will expand (or shrink) exponentially in the 
grand-canonical simulation. By requiring an exponential increase (or decrease) to 
disappear in the simulation we can find the parameter A precisely. We then generate 
configurations dynamically near to the critical point Ac, constraining N2 within a 
upper- and a lower-bound (A^™" < N2 < N^°-^). 



2.2 Procedure of the grand-canonical simulation 

We generate dynamically triangulated surfaces using the grand-canonical Monte 
Carlo method. In order to keep N2 within the appropriate range wide enough to 
extract 7 from eq.(|]), but not too large for the sake of the computing time, we 
simply ignore all of the moves that make the number of triangles become less than 
]\fmzn greater than N^"-^. Here, we restrict ourself to graphs excluding "tadpole" 
and "self-energy" in the dual lattice. We note that in the grand-canonical method 
we must pay careful attention to satisfying the detailed valance, because the number 
of states depends on the total number of triangles and the coordination number of 
each vertex. 

We first explain the method used to generate dynamically triangulated surfaces 
in the case of pure gravity (c = 0). Suppose that we wish to increase the number 
of real vertices from Nq to Nq + 1; to do this we choose two dual links belonging 
to the same dual loop, and insert a "propagator" dual link joining the two dual 
links, thus creating an additional dual loop. For the inverse move, Nq Nq — 1, we 
choose a dual link at random, and remove it, resulting in two dual loops separated 
by the dual link merging into one dual loop. The total number of possible increase 
(or decrease) moves from configuration A are 

i=i ^ 

No „. 

nT = T.% (11) 

i=l 
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Table 1: Sphere 



h= 


=0 


7 


Ac 


Pc 


c= 




-0.517(27) 


1.12467 




c= 


:0.5 


-0.345(27) 


1.53716 


0.227 


c= 


:1.0 


-0.009(26) 


1.94848 


0.226 


c= 


:1.5 


0.112(27) 


2.35843 


0.225 


c= 


-2.5 


0.238(26) 


3.17513 


0.223 



where {qi} is the set of coordination numbers of configuration A. 

In these moves the ergodic property is automatically satisfied. However, we must 
treat the detailed valance carefully, since the total number of possible moves depends 
on Nq and {qi}- We must multiply the prefactor to the Boltzmann weight in order 
to satisfy the detailed valance equation for the moves between configurations A and 
B, given by 

-^^P (-12) 

where Pa^b is the transition probability from configurations A to 5, Wa is the 
Boltzmann factor, and ha is the number of possible moves from configuration A. 
For detailed arguments for eq. ([T^ ) on random lattices, see ref . [jll] . 



For the case d ^ when a dual loop is separated by two dual loops, two spins 
should be attached to the two dual loops, which counts four spin states; in the 
inverse case one of two neighboring spins should be taken out, which counts two 
spin states. The combined number of operations and the Boltzmann factor are 



nA = y^{ ' ■ 4'^ + ^ ■ 2^^ , (13) 




Wa = exp { -XN2 + /? E E \ ■ (14) 

M=l <ij> 



For the spin variables, Monte Carlo trials with Wolff's cluster algorithm ||T2| are 
performed after some geometrical moves. 



3 Numerical results 

We have measured the string susceptibility exponent for various combinations of 
c = I = 0, |, 1, |, I and h = 0,1, 2. The central charge and topology of the surface 
are kept fixed during each simulation. When Ising spins are put on the surface, we 
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Table 2: Torus 



h= 


=1 


7 


Ac 




c= 




2.014(24) 


1.12469 




c= 


:0.5 


1.995(23) 


1.52890 


0.2163 


c= 


:1.0 


1.971(22) 


1.93335 


0.2163 


c= 


:1.5 


1.980(23) 


2.33841 


0.2163 


c= 


-2.5 


2.018(25) 


3.14950 


0.2163 



Table 3: Double — torus 



h= 


=2 


7 


Ac 


Pc 


c= 


:0 


4.510(27) 


1.12471 




c= 


:0.5 


4.338(26) 


1.52821 


0.215 


c= 


:1.0 


4.043(27) 


1.93039 


0.214 


c= 


:1.5 


3.909(26) 


2.33123 


0.213 


c= 


:2.5 


3.781(27) 


3.12972 


0.211 



search the critical coupling strength (/5c) with the fixed number of triangles set at 
3000 before making a grand-canonical simulation. We determine the critical coupling 
strength from the peak of magnetic susceptibility for each combination of c and h. 
In the case of a torus the experimental values (/3c) were quite close to the theoretical 
value (Pc ~ 0.2163), and we used the exact value; in the other cases, however, we 
used the experimental Pc- We set the maximum number of triangles [N^"-^) to 3000 
and the minimum one (iV™") to 1000. One sweep is defined by 2000 x 2'^ trial moves 
and 5 cluster operations of Ising spin. 

In order to obtain P{N2), we measure the number of triangles (A^2) for 10^ times 
every 5 sweeps. Plotting lnP(A''2) versus lnA^2; such as Fig.[T|, we determine the 
string susceptibility exponent from its slope. Our results are summarized in Tables 
1, 2 and 3 for a sphere, torus and double-torus, respectively. Especially, in the pure 
gravity case we know the theoretical value, 

1 256 

Ac = -In— ^ 1.12467, (15) 

and our experimental Ac is very close to theoretical prediction. The errors indicated 
in parentheses for the last two digits are only statistical errors estimated from the 
least-squares fits. The systematic errors are due mainly to how close we can de- 
termine A to Ac. In our simulation this systematic error is estimated to be around 
0.02 ~ 0.03. 
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Figure 1: Distribution (P(A^2)) for spherical pure gravity with log- log scales. The 
parameter A is 1.12467, and the string susceptibihty exponent (7) is —0.517. 
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Figure 2: String susceptibility exponents (7) of a sphere, torus and double-torus 
for pure gravity {c — 0). 
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Figure 3: String susceptibility exponents (7) of a sphere with various central charges 
(c). Solid line represents eq.(2) with h = 0. 
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Figure 4: String susceptibility exponents (7) of a torus with various central charges 
(c). Solid line represents eq.(2) with h = 1. 
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Figure 5: String susceptibility exponents (7) of a double-torus with various central 
charges (c). Solid line represents eq.(2) with h = 2. 

4 Summary and discussion 

We calculated the string susceptibility exponents of 2-dimensional triangulated sur- 
faces coupled to d-Ising spins ((i=0,l,2,3,5) with various topologies (/i=0,l,2) using 
the grand-canonical Monte Carlo method. One of the advantages of our method is 
the ability to treat various topologies. Fig.^ shows good agreements to the theory 
for not only a sphere, but also for a torus and a double-torus of the pure grav- 
ity. Furthermore, in Figs.^^ and ^ the string susceptibility exponents for various 
topologies {h = 0, 1, 2) are reproduced very closely to the theoretical predictions as 
long as c < 1. Especially for c = 1 of the spherical surface we give the theoretically 
expected 7. The log correction in the partition function is one order smaller than 
the leading term in our simulation for N2, being on the order 10^, which proves that 
our method is free from a finite-size correction, compared to the MINBU algorithm. 

We are especially interested in the behavior of 7 for c > 1. In the case of a torus 
{h = 1) we cannot observe the effect of the central charge in Fig.|, as expected 
from the theoretical prediction in eq.(^ and eq.(|^). In the spherical case {h = 0) 
7 becomes positive and increases slowly, while in the double-torus {h = 2) it seems 
to decrease slowly. Although 2-dimensional surfaces seem to make a transition 
to the branched polymer phase gradually as c > 1, based on these figures, we 
cannot conclude for certain what kind of branched polymers were obtained in our 
simulations. It is preferable to find other quantities to observe the branched polymer 
phase more clearly. For example, in the case of a higher genus the moduli of the 
surface is one of the most important quantities in order to distinguish the branched 
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polymer |jT3 . 
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